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Abstract 

We analyze the general nonlinear evolution equations for multi gluon correlators derived 
in |jj by restricting ourselves to a double logarithmic region. In this region our evolution 
equation becomes local in transverse momentum space and amenable to simple analysis. 
It provides a complete nonlinear generalization of the GLR equation. We find that 
the full double log evolution at high density becomes strikingly different from its linear 
doubly logarithmic DGLAP counterpart. An effective mass is induced by the nonlinear 
corrections which at high densities slows down the evolution considerably. We show that 
at small values of impact parameter the gluonic density grows as a logarithm of energy. 
At higher values of impact parameter the growth is faster, since the density of gluons is 
lower and nonlinearities are less important. 
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1 Introduction 



One of the most important theoretical problems in the high energy hadronic physics is the un- 
derstanding of high density nonlinear effects which may lead to unitarization of perturbatively 
generated growth of the hadronic cross sections at very high energy. The simplest process in 
which one hopes to see these effects is deep inelastic scattering. Standard linear evolution Q], 
which is phenomcnologically so successful at present energies predicts a perturbatively generated 
growth of the total cross section faster than a power of a logarithm. This contradicts the expecta- 
tion based on the unitarity constraint that the growth should not be faster than the second power 
of the logarithm of energy^], ft is natural to assume that the effects that are left out of the linear 
evolution are the ones that are ultimately responsible for the unitarization effects. In particular 
the suspicion falls on possible screening effects due to finite partonic density. At high energy this 
density grows and the emission of high transverse momentum partons which are scattered by the 
probe should be inhibited relative to the low density situation []. These nonlinear effects are not 
taken into account by the linear evolution equations. Recently it has been argued that already 
the presently available data on DIS can not be reasonably explained without taking into account 
the nonlinear effects in QCD evolution ||. 

Possible nonlinear generalizations of perturbative evolution have been discussed many times 
in literature starting with the famous GLR model || . A fully nonlinear evolution equation based 
on a particular model for the gluon source was derived by Mueller in ]10[ |, which inspired much 
of the later research in this field. More recently generalized nonlinear evolution equations were 
discussed by Levin and Laenen and Ayala, Gay-Ducati and Levin |L2| . All these approaches 
are based to some extent on a physically motivated ansatz and it is important to understand the 
impact of nonlinearities in a more general framework. 

In a recent series of papers p|, p^I5| we have developed an approach to the evolution of 
dense partonic systems within the framework of the Wilson renormalization group^J The method 
is based on earlier work by McLerran, Venugopalan and Ayala, Jalilian-Marian, McLerran and 
Venugopalan p9| in which the idea of representing fast partons by static color charge density was 
developed. This approach results in a nonlinear functional evolution equation for the generating 
functional of the color charge density correlators, which is valid to leading order in a s at densities 
which parametrically do not exceed l/a s . This is indeed the region of density which is interesting 
for the perturbative unitarization effects. The equation is fairly complicated since it only requires 
ordering in longitudinal momenta during evolution and puts no constraint of any kind on the 
ordering of transverse momenta. In fact, in the low density limit it reduces to the BFKL equation 
for the two point correlation function fl4|| and a simple truncation of it reproduces the complete 
BKP hierarchy |(J. 

The aim of the present paper is twofold. First we are going to show how by a simple trans- 
formation this evolution equation is converted into a set of equations for the evolution of the 
correlators of the chromoelectric field rather than the color charge density. This form is perhaps 
more advantageous, since these correlators are more directly related to observable quantities and 
for example the DGLAP evolution operates directly with them. Our second goal is to consider 
these evolution equations in the doubly logarithmic regime. That is to say we will work in the 
approximation where strict transverse momentum ordering is imposed on the evolution. In our 

1 There is a certain caveat here, since in a strict sense there is no proof that the DIS cross section has to unitarize 
in the same way as purely hadronic cross sections. Still one believes based on pjwsical arguments that the growth 
of the DIS cross section at high Q 2 should not exceed a power of a logarithm [H, 14] - 

2 Another possibility is that even at high Q 2 at high enough energy the scattering becomes nonperturbative 
due to large contributions from the small transverse momentum region. This situation then is entirely outside the 
reach of the methods of perturbative QCD. It has been however convincingly argued that even if this is the case 
asymptotically, there are many physically interesting situations where the main contribution to the cross section 
come from the perturbative region Q-Q]. In these cases the perturbative nonlinearities must play a crucial role. 

3 Let us note that there exists a number of approaches to small-x physics based either directly on constructing 
the QCD effective action at small x by combining the reggeon and usual gluon terms |16| ] (see also Wm), or on the 
generalization of the operator product expansion [ |lq ] . It is important to understandthe relation between these 
different approaches and the approach described in this paper. At the moment this relation is not completely clear 
to us. 
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framework this corresponds to the leading order in the expansion in powers of transverse deriva- 
tives. It turns out that the evolution equations simplify tremendously in this limit and become 
much more tractable. 

The main feature of the full nonlinear evolution is the appearance of a dynamical mass. This 
mass is induced by quantum corrections to the QCD evolution at finite density and is proportional 
to the square of the chromoelectric field. Its effect is to slow down the evolution relatively to the 
standard perturbative case. As a consequence the gluon density in the small impact parameter 
region at high density grows with energy only logarithmically as opposed to the exp Vhos type 
of growth at low density. To construct an explicit solvable model we approximate the generating 
functional for the chromoelectric field by a Gaussian. In this Gaussian approximation the only 
independent quantity is the distribution function itself, that is the two point correlation function 
of the chromoelectric field. All higher correlators have simple factorization properties in terms of 
the two point function. This model leads to a simple closed nonlinear equation for the distribution 
function which is qualitatively quite similar to the equation considered in fl2]| . 

The paper is structured as follows. In Sec. 2 we briefly review the general structure of the 
nonlinear evolution equation and show how to rewrite it directly in terms of correlators of the 
chromoelectric field. In Sec. 3 we explain how the double logarithmic limit arises in the present 
framework and derive the explicit evolution equation in this limit. Sec. 4 is devoted to analysis 
of the the resulting equation and the derivation of the simple Gaussian model. Finally Sec. 5 
contains a brief discussion. 



2 Evolution equations for the correlators of the chromo- 
electric field. 

First let us briefly recall the framework and the results of Q. In this approach the averages of 
gluonic observables in a hadron are calculated via the following path integral 



<0[A] >= J DpDA„0[A]exp(- J d 2 x ± F[p a 



(s-0] (1) 



d^x-F^F^-— j d 2 x ± dx~5(x-)p a (x ± )trT a W^ 00 , 00 [A-}(x-,x ± ) 



ipFp" L 

'4 a N c 

where W is the Wilson line in the adjoint representation along the x + axis. The hadron is 
represented by an ensemble of color charges localized in the plane x~ = with the (integrated 
across x~) color charge density p(x±). The statistical weight of a configuration p(x±) is 

Z = cxp{-F[p}} (2) 

In the tree level approximation (in the light cone gauge A + = 0) the chromoelectric field is 
determined by the color charge density through the equations 

F +i = U{x-)a i {x A _) (3) 

and the two dimensional vector potential on (x ± ) is " pure gauge" and is related to the color charge 
density by 

did} - - r hc a\a) = 

= ~p a (4) 

Integrating out the high longitudinal momentum modes of the vector potential generates the 
renormalization group equation, which has the form of the evolution equation for the statistical 
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weight Z §[T| 



2 5p(u)Sp(v] 



[Zx(u, v)] 



8p(u 



[Za(u)] 



(5) 



In the compact notation used in Eq. (||), both u and i> stand for pairs of color index and transverse 
coordinate, with summation and integration over repeated occurrences implied. The evolution in 
this equation is with respect to the rapidity £, related to the Feynman x by 



C = ml/a; 



(G) 



Technically it arises as a variation of Z with the cutoff imposed on the longitudinal momentum 
of the fields A^. The quantities x[p] an d cr[p] have the meaning of the mean fluctuation and the 
average value of the extra charge density induced by the high longitudinal momentum modes of 
A^. They are functionals of the external charge density p. The explicit expressions have been 
given in |j| . For the purpose of the general discussion in this section we do not need their explicit 
form. Eq. (|^) can be written directly as evolution equation for the correlators of the charge density. 
Multiplying Eq. (|^) by p(xi)...p(x n ) and integrating over p yields 



<K 



< p{xi)...p(x n ) > 



(7) 



a s 



E 

.0<m<fe<n+l 



E 

0<i<n+l 



< p(x 1 )...p(x m -i)p{x m+1 )...p(x k - 1 )p(x k+ i)...p(x n )x(x m ,Xk) > 
< p(xi)--p(xi^ 1 )p(xi +1 )...p(x n )a(xi) > 



In particular, taking n = 2 we obtain the evolution equation for the two point function 
— < p{x)p{y) >= a s {< x(x, y) + p(x)a(y) + p(y)cr(x) >} 



(8) 



This set of equations for the correlators of the color charge density completely specifies the 
evolution of the hadronic ensemble as one moves to higher energies (or lower values of x). 

Our aim in this section is to convert Eq. (^|) into a set of equations which determine directly 
the evolution of the multi particle correlators of the chromoelectric field ct; ^ For this purpose we 
multiply Eq. (||) not by p{xi)...p{x n ) but rather by a" 1 (xi)...a^ n (x n ) and again integrate over p. 
This results in a set of equations analogous to Eq. (||) 



<< i 1 (x 1 )...< : (x„)> 



0) 



0<i<n+l 

E 

0<m<fc<n+l 



< <(x 1 )...<: i i (^_ 1 )< ! + + l i (x ;+1 )...< : (x„)< ! (^) > + 



Here we have defined 



a k-l / \ a fc + l / \ flu/ \ a-rndk I \ ~~ 

xt 3 b ( Xl y) >= r? c (x, u) X cd (u, v)rf\v, y) 
<(x) = rf(x, u)a b (u)+pf c (x, u, v) X hc (u, v) 



(10) 



4 A11 the functions in the rest of this paper depend only on transverse coordinates. To simplify notation we 
largely drop the subscript _L in the following. 

5 The chromoelectric field strength F 1+ is proportional to a; as per Eq. (^). In what follows we will therefore 
refer to on as to the field strength. 
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with 



r?(x,y) 



5p b (y) 



(11) 



pf c (x,y,z) 



5 2 a?{x) 
6p b {y)6 P c (z) 



In Eqs. ([Tl]), ( |l2] ) the repeated indices are summed over. 

We now have to express the quantities r and p explicitly in terms of the field a, . Once this is 
done the equations Eq. ([lO]) become explicit equations for the evolution of the correlators of the 
field cti, since the charge density p is already known in terms of on by virtue of Eq. (^). 

To find rf b we differentiate Eq. (0) with respect to p and get 



dtrf = -5 ab 



(12) 







This set of equations is easily solved by decomposing into the longitudinal and the transverse 
parts according to 

rf = d,l ab + e lJ d j t ab 



(13) 



with the following result 

rf(x,y) = -<x\{d i 6 ab + e ij d j 



Dd 



£kia c k di } -^\y >= - < x\ 



-i ab 



D, 



dD 



\y > (14) 



Here < a;|0|y > denotes a configuration space matrix element in the usual sense. Scalar products 
with respect to space time indices here and below refer only to transverse indices, so that for 
example dD — diDi with i = 1,2. For convenience we have also defined 



af 



fab c a c 

d z 6 ab + af . 



(15) 



Now differentiating Eq. ( |l3| ) once again with respect to p we can solve for p: 

1 



pf c (x,y,z) 



eijdj 



Dd 



ad\ 



(x,u)f d f e e kl r f k b (u,y)rr(u,z). 



(16) 



Here again summation over the repeated color indices and integration over the transverse coordi- 
nate u on the right hand side is understood. 

Equations ([u]) and ([If]) give the complete solution for r and p in terms of the field a. The set 
of equations Eq. ( |l0| ) together with Eq. ( pd| ) and Eqs. (|l4|), ( |l6| ) directly govern the evolution of 
the correlators of chromoelectric field. 

We close this section by noting that the quantities xtj an( i a i have very simple physical 
meaning. As mentioned above, the high momentum modes of the vector field which have been 
integrated out in order to arrive at the evolution equation induce extra color charge density 5 p. The 
average value of this induced density and its mean fluctuation appear in the evolution equations 
for the correlators of charge density, as a a and x ab ■ Clearly the appearance of the induced color 
charge density leads to the change in the value of the chromoelectric field through the solution of 
Eq. (^) with p + 5p on the right hand side. Diagrammatically the new field is represented by the 
sum of the tree level diagrams with p + Sp as the source. Those are depicted in Fig. [l[ It is clear 
that <rf is nothing but the average value of the induced chromoelectric field. Additionally, since 
the induced charge density Sp contains components with frequencies p~ higher than those of the 
input charge density p, upon time averaging the induced chromoelectric field is characterized also 
by mean fluctuation < 5afSa b >= xtj ■ This again has the simple diagrammatic representation 
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Figure 1: The diagrams contributing to the chromoelectric field at order a s . The full circles denote the 
background charge density p. The empty circles denote the average of the charge density induced by the 
fluctuations < 5p >= a. The black bars denote contractions corresponding to the mean fluctuation of the 
induced charge density < SpSp >= X- 




Figure 2: The average fluctuation of the induced chromoelectric field xtj ■ The notations are the same 
as in Fig.l. 

shown in Fig |2|.[] The explicit expressions Eqs. ([n]), (fl4| and ( [l6| ) are equivalent to resumming all 
the diagrams of the type indicated in Fig. |] and Fig. |2j. They clearly exhibit the mechanism of 
the evolution of the field correlators in our approach. 

3 The Doubly Logarithmic Regime 

We would now like to consider the evolution in the doubly logarithmic approximation. By this 
we mean that in every step of the evolution not only the longitudinal momenta are lowered but 
in addition, the transverse momenta are required to grow. That is of course just the QCD parton 
model picture that the "valence" partons have transverse momenta on the typical hadronic scale 
while the high transverse momentum partonic components of a hadron exist only as fluctuations 
at short time scales. As discussed in jjj increasing the frequency in the evolution is equivalent to 
transverse momentum ordering in addition to the longitudinal momentum ordering. This scheme 
is therefore equivalent to the Born-Oppenhcimer approximation in a system with a continuum of 
time scales. 

Technically, in our approach the doubly logarithmic regime is singled out by restricting the 
transverse momenta in the fluctuation field which is being integrated over to be smaller than the 
transverse momentum of the observable that is being calculated, but larger than the transverse 
momentum of the background field. Practically this just means that the background field through- 
out the calculation is considered to be independent of x±. The results are then easily obtained 
by taking the leading order of the transverse derivative expansion of the general expressions for \ 
and a given in [jjj. The result in this limit simplifies very much for two main reasons. First, since 
there is ordering in all momentum components in the evolution it is clear that virtual diagrams 
can not give a non-vanishing contribution. This is simply the result of momentum conservation 

"Given that Eq. (Jsjj contains the average value and the mean fluctuation turn as the only relevant characteristics 
of the induced density 8p this naturally carries over to the the induced field 5a to leading order in the coupling 
constant a s . 
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in the xj_ independent background. Therefore we have immediately 

a = a, ■, = U . 



(17) 



The calculation of the real part also simplifies. This is because the two components of the 
chromoclectric field commute in this limit 



[cm, atj] = 



(18) 



This is the consequence of the first equation in Eq. (|4|) for constant cm. This means that cm can be 
treated as constant numbers. Using the explicit expressions from 0] for x we obtain in the double 
logarithmic limit 

(3D) 2 



and 



a < x\ 



D 2 (d 2 + 2a 2 ) 



\y > 



a < x\ 



DiDj 



■\y > 



ab 



(19) 



(20) 



, D 2 (d 2 + 2a 2 ) 

Let us first check that our expressions in the limit of weak field do indeed reproduce the 
DGLAP evolution in the double logarithmic regime, as claimed. To this end let us consider the 
evolution of the two point function^ . 



— < <(x)a°(y) >=4a s 



< x\a 



DiDi 



D 2 {d 2 + 2a 2 



■\y > 



(21) 



At weak fields the 2a 2 term in the denominator of Eq. (^l|) should be dropped and the covariant 
derivatives become simple derivatives. Taking trace of the evolution equation over both, color and 
transverse indices we obtain 



-| < at{x)a a i {y) >= -4a s 7V c < « >< x\^\y > 
The gluon distribution is related to the field correlator in the following way Q| 
xG(Q) = (2n f J d 2 b< F +l (b,xp + ,x+)F +l (b,xp + ,x + ) > Q 



(22) 



d 2 b 



d z k I d 2 xe lkx < a?(b + x/2)a?(b - x/2) > 



(23) 



Here F +l (b,k + ) — f dx e lk+x F +l (b,x ) and the subscript Q means that the transverse coordi- 
nates in the chromoelectric fields operators are equal with accuracy 1/Q. This also gives 







d\nQ- 



\xG(Q)] = -^-Q 2 J d 2 b J d 2 xe lQx < a?(6 + x/2)o?(6 - x/2) > 



(24) 



The factor — appears in these expressions since the vector potential on differs in normalization 
from the chromoelectric field by a factor 1/g, see Eq. (||). Taking the Fourier transform of Eq. (|2j 
and integrating over the impact parameter bi we obtain 



(fi 



dQdlnQ- 



[xG{Q)\ = ^^ X G(Q) 



(25) 



which is precisely the double logarithmic approximation to the DGLAP equation. 

7 In this equation as well as in the rest of the paper we rise the bra and ket notation < | > to denote both the 
averaging over the hadron state (averaging over p in Eq. An)) of the function of a^'s and the configuration space 
matrix element of the powers of derivatives di . 
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The following comment is due here. The right hand side of Eq. (|22|) contains the expression 
< afa" >. In the leading order of the derivative expansion on does not depend on the transverse 
coordinate, and this expression of course is also constant. If one allows slow variation of the field 
(slow in the sense that the scale of the spatial variation should be much larger than 1/Q) this 
obviously should be understood as < af(b)af(b) > where the impact parameter b = (x + y)/2 
is the center of mass coordinate. Of course since the spatial resolution is 1/Q, the coordinate of 
the field is only defined to this accuracy. The more exact specification can be only given beyond 
the leading order of the derivative expansion and is therefore beyond the approximation we are 
discussing here. Another way of saying this is that the fields that contribute to the right hand 
side of the evolution equation themselves contain only transverse momentum components smaller 
than Q. We therefore have to understand < a°a° > on the right hand side of Eq. (p2|) as 

< <(&)<(&) >= £ 0| d 2 xe lkx < a«{b x/2)a?(b + x/2) > (26) 

This is precisely how G(Q) arises in the right hand side of Eq. (p5j). In this approximation therefore 
the dependence on Q and b in Eq. (^2|) factorizes and the products of the fields are defined with 
the transverse cutoff 1/Q- This is a general feature of the leading order of derivative expansion 
and is not specific in any sense to the weak field limit. 

We note that expanding the right hand side of Eq. ( ^l[ ) to second order in the field intensity 
and taking trace over color and Lorentz indices we obtain the first nonlinear correction to the 
DGLAP equation 

ih - xG {Q) = ^xG(Q) - A f d 2 bti < (a(b) 2 ) 2 > (27) 



<9C<91nQ 2 TT Q 2 

This is precisely the contribution of the twist four operator to the evolution of gluon density as 
calculated by Mueller and Qiu ^J. It * s probably worthwhile noting that Eq. ( p7| ) should be 
compared not with the nonlinear GLR equation, which is expressed in terms of the distribution 
function, but rather with Eq. (25) of plj which contains the contribution of twist four operators to 
the evolution of the gluon distribution. It is straightforward to check that taking into account the 
normalization of the operator G' 2 ' in pl[ our Eq. ( p7| ) is indeed equivalent to Eq. (25) of plfl . The 
nonlinear GLR equation arises as an approximation to this equation if one assumes factorization 
of the four gluon operator into the square of the gluon distribution in a very particular manner. 
Assuming this factorization Eq. ( p7j ) would indeed lead to the GLR equation with the same 
coefficient as given in pl| . 

Equations for the higher correlators can be analyzed in a similar fashion in the weak field 
limit. They all reduce to homogeneous evolution equations and yield simple expressions for the 
anomalous dimensions of the higher twist operators. Our main interest here, however, is to explore 
the opposite situation, namely when the fields and partonic densities become large and this is what 
we will do in the rest of this paper. 

The most prominent feature of equation (|2l|) is the appearance of the factor Q 2 — 2a 2 in the 
denominator. The field strength —2a 2 therefore plays the role of a dynamically generated mass 
in the evolution at high density^. This is not entirely unexpected, but it is gratifying to see the 
emergence of this dynamical mass in such a straightforward fashion in our approach. Clearly the 
effect of this mass is to slow down the evolution. In fact it is easy to see that when the density 
becomes large this "mass term" dominates the evolution and leads to unitarization. 

Let us consider the evolution of the distribution function in the high density -strong field limit 
—a 2 > Q 2 . Neglecting Q 2 relative to a 2 in the denominator in Eq. (21)^] we obtain 



° 2 ntn\ - N c(N c - 



dCdl,Q 2 ^ Q) = ^^ QS (28) 

8 Note that the matrix a as we have defined it is anti hermitian. The eigenvalues of a 2 are therefore neg ative. 



This must be done with care, since the matrix a has zero eigenvalues. See comment following Eq. (|2q). 
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Here 5* is the area of the hadron within which the density is high. It appears due to the integration 
over the impact parameter. The origin of the color factor N C (N C — 1) is the following. If the 
matrix a 2 had no zero eigenvalues, this factor would be just N 2 — 1 - the trace of the unit matrix 
in the adjoint representation. However the matrices on as defined in Eq. ( |l5[ ) necessarily have zero 
eigenvalues. The number of these zero eigenvalues is generically N c — 1. Since [a,*, otj] — 0, the same 
applies to the matrix a 2 . The contribution to the trace therefore comes only from the subspace 
spanned by the eigenvectors which correspond to nonzero eigenvalues. The dimensionality of this 
subspace is N 2 - 1 -JN C - 1) = N C (N C - 1). 



Integrating Eq. (|28|) we obtain 

X G(Q 2 ) = iVc(A :~ 1) Q 2 ffhi(l/.T) + Go (29) 



This result is very interesting. First it shows a characteristic feature expected of unitarization - 
the growth with energy is only logarithmic rather than power like as in the case of a pomeron, 
or exp ^/\og s as in the case of the doubly logarithmic DGLAP evolution. On the other hand 
the gluon density does not completely saturate but keeps growing logarithmically even at high 
densities. Similar "partial saturation" was found in the analysis of p0| . and also in the solution 
of a nonlinear evolution equation in jl2j . 

Identifying naively the gluon distribution with the total DIS cross section for a probe that 
couples directly to gluons^ one would get a logarithmic growth of the cross section on top of 
the factor of the geometric area of the hadron. This is perfectly compatible with the Froissart 
bound and in fact grows slower than the allowed square of the logarithm if the effective area S in 
Eq. ( p9| ) grows slower than logarithmically. On the other hand the prefactor depends strongly on 
the number of colors and is large at large N c . This is not what one expects generically from the 
unitarized cross section. Once the geometric size has been factored out one does not expect any 
other numerical factors. The factor N 2 is of a purely perturbative origin. It counts the number of 
perturbative degrees of freedom - gluons that take part in the scattering. This is in accord with the 
fact that our approach is indeed perturbative and we therefore can only detect the perturbative 
mechanism for unitarization which should be operative in the region where the scattering is mainly 
perturbative. Apart from the perturbative component which is described by our present approach 
there is also a nonperturbative, low transverse momentum one of the soft pomeron type (see e.g. 
the recent analysis in p2[|) which is not governed by the evolution we consider here. At very low 



values of x it must become leading and eventually the perturbative behavior of Eqs. (|28[), (|29|) will 
cross over into the real asymptotics which may behave as a square of the logarithm. 



4 The Gaussian model. 

Eq. (|2^) represents the asymptotics of the distribution at large values of the density. At small x 
one expects the density to be large close to the center of a hadron, at small values of the impact 
parameter, but not near the boundary. Eq. ( |28| ) therefore can not be quite correct even at small 
x. To get a more realistic picture one has to allow for a possibility of a more rapid evolution in 
the peripheral regions. 

Generically the system of equations Eq. ([h]) is a coupled system of equations for infinite 
number of the multi gluon correlation functions and as such is nontrivial even in the doubly 
logarithmic limit. We would like to simplify this system so that it becomes more tractable. One 
straightforward possibility is to consider a Gaussian truncation of the original system. What we 
mean by this is to assume standard Gaussian factorization of multi particle correlators of the type 

< a 2n >~ [2n- 1)!! < a 2 >" (30) 

10 This naive identification has to be taken with a grain of salt. For large values of the gluon density Eq. ( ps| ) the 
leading twist relation between the cross section and the number of gluons is not valid. One then has to consider the 
contribution of higher twist operators directly in the cross section as well as in the evolution of the gluon density. 
The following argument should therefore be understood only in illustrative sense. 
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This approximation is in the spirit of the GLR model || . 

Technically this means that we take the statistical weight used for calculation of the averages 
on the right hand side of the evolution equation as a Gaussian function of the fields. One should 
of course take care to ensure that this Gaussian weight is consistent with the symmetries of the 
system and also with the fact that the two spatial components of the chromoelectric field are not 
completely independent but rather commute with each other, see Eq. (|l8|). To ensure this, let 
us explicitly solve the constraint commutativity imposes on the matrices on. The quantities to 
be averaged all being invariant under global color rotations implies they can only depend on the 
eigenvalues of the on. Since the on commute, they can be diagonalized simultaneously. We can 
therefore take the following explicit representation^ 

ai =iJ2t?T A (31) 

A 

Here T A are the diagonal Cartan subalgebra generators of the SU(N) group in the adjoint repre- 
sentation, tf are real coefficients and the index A takes values from 1 to N c — 1. The coefficients 
tf are completely independent and the Gaussian weight for the averaging can be taken as 

< ... >= nj^)-' [n„ Adt ? ... e*p{-fc±2 f ^} (32) 
V a s xg J J I a s xg ) 

The meaning of the parameter g is made clear by calculating the glue-glue correlator 

<«?<>=%• (33) 

The quantity g is therefore the local partonic density. Just as the field itself it is of course a slowly 
varying function of the impact parameter, and is related to the gluon distribution by 

G(x, Q 2 ) = J d 2 b ± g(x, b ± , Q 2 ) (34) 

Fourier transforming Eq. ( |2l| ) and taking trace we get for the right hand side 

Ar - _1 fAj-B A B 

where rf is the n-th eigenvalue of the A-th Cartan subalgebra generator in the adjoint represen- 
tation. The N c — 1 dimensional vectors r„ are the root vectors of SU(N) and have the following 
important properties 

v 2 n = 0, n = l,...,N c -l (36) 
r 2 n = 1, n = N Cl ...,N 2 c -l 

Therefore only the terms n = N c , ...,N% — 1 contribute to the sum. In each one of these non 
vanishing terms one can perform an orthogonal rotation of the coefficients tf — > sf such that 

si = tfrf (37) 

The possibility of such an orthogonal rotation is assured by the fact that the root vectors are 
properly normalized. The expression in Eq. (p5| ) then simplifies to 

dsi exp \ s i f n.i I n„2 (38) 

2 



xg J { oi s xg 1 ) Q 2 + 2s 



nN c (N c - 1) 2 Q 2 



xg Jo 



d 



s 




l + s_ 


exp | 





AN C - 1)Q 
2a s xg 



11 The most general matrix a; has the form iU ^] .t^T^U^ , with U an adjoint representation SU(N) matrix. 
However, as noted before any SU(N) invariant quantity does not depend on U and we therefore omit it in all the 
following formulae. 
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The evolution equation then becomes 



d 2 

xg(x,Q,b±) (39) 



d(dlnQ 2 

N C (N C - 1 



Q~ 



. ^(7V C -1)Q 2 /7r(iV c -l)Q 2 \ / ^V C -1)Q 2 
1 H exp < > ca 



2a s xg \ 2a s xg J \ 2a s xg 

where Ei(x) is the integral exponential function defined as 



Ei(a;) = - / dt , x < (40) 

•J — x t 

We stress that the factorization assumed in Eq. (|32) is motivated only to the extent that it 
leads to a closed evolution equation for the gluon density. One can think of it as a mean field 
approximation for the true weight function for the averaging over fields a, . One certainly expects 
this approximation to be valid at small xg where it is equivalent to the steepest descent integration 
over afj's. Additionally at large xg the asymptotics of Eq. ([59]) reduces to the exact asymptotics 
of the original equation (|2l|) and is independent of the form of the weight function. It is therefore 
reasonable to expect that this simple model provides a sensible interpolation of the evolution 
between the weak and strong density regimes. 



Note that Eq. (|39| ) governs the evolution of the gluon density xg(b) and is local in the impact 
parameter space. It therefore contains more information on the structure of the hadron state 
than just the gluon distribution xG(x, Q 2 )- To determine the evolution of the gluon distribution 
xG(x,Q 2 ) one has to integrate the solution of Eq. ( ^9|) over the impact parameter b. To write 
down a closed equation directly for xG one would have once again to resort to modeling. This 
time one needs a model of the gluon density as a function of the impact parameter. Models of 
this kind have indeed been used in the literature. One usually assumes a simple factorizeable 
structure Q g(x,Q 2 ,b) = G(x 1 Q 2 )S~ 1 (b). The resulting equation for G then depends on the 
form of 5(6) one takes in this ansatz. For example for a simple homogeneous disk model 5 _1 (6) = 

{irR^^eiR 2 - b 2 ), Eq. Q with the substitution Q 2 Q 2 ttR 2 governs the evolution of xG. 

_ b 2 

More commonly one uses the Gaussian type distribution 5 _1 (6) = (irR 2 )~ 1 e 'zr 1 . Any specific 
ansatz of this type introduces a new phenomenological parameter R into the equation which has 
the meaning of an effective radius of the gluon distribution in the hadron. The same kind of 
parameter appears in the GLR equation 0, fH|. 

It would be interesting to study equation ( |39|) numerically to find out how the effective area in 
which the partonic density is high [S in Eqs. (|28|), (p9[)) depends on the rapidity. This we intend 
to do in a future publication. 



5 Discussion 

We have considered the doubly logarithmic limit of the low x nonlinear evolution equations derived 
in jij. The most salient feature in this regime is the appearance of a dynamically induced mass 
which leads to the slowdown of the evolution as the partonic density becomes larger. In particular 
we have shown that at small impact parameters the gluon density grows as a logarithm of energy 
rather than following the much faster growth pattern predicted by DGLAP evolution. To study 
the crossover between these two regimes in detail one has to analyze the set of equations ( |l0| ) for 
the evolution of multi gluon correlators. This is still a formidable problem, although in the doubly 
logarithmic limit the equations are much simpler than the general evolution equations of |l| . 

We have discussed an approximation in which a closed equation Eq. (|3|) for the gluon density 
arises. This approximation imposes Gaussian factorization of multi particle operators in terms 
of the gluon density and is similar in spirit to the approximation used to write down the GLR 
equation in and pi]]. 



We note that the local equation for the gluon density Eq. (39) is qualitatively similar to the 
equation suggested in Jl^| . The authors of [|12| suggested the following equation on the basis of 
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reinterpretation of the Glauber formula for multiple scattering of |K 



if 



d(d\nQ' 



-xg 



2Q 2 



(41) 



Although the functional form of this equation looks somewhat different from our Eq. ( p^ ) the 
main features are the sameQ. Both reduce to the DGLAP equation in the limit of small density. 
Asymptotically at low x, where the density is high the solution to both equations behaves as 
Q 2 lnl/x. 

There are however also important differences between the two equations. At large densities 
the correction to the asymptotic solution of Eq. ( p| ) is exponentially small. On the other hand 
expanding Eq. (p9|) to order 1/xg we find 



dCdlnQ' 



: xg(x,Q,b) 



N C (N C - 1) Q2 



l ^V c -l)Q 2 ln 
2a s xg 



2a s xg 



n(N c - 1)Q 2 



(42) 



So the correction here is a power enhanced by a logarithmic factor. The approach to the asymp- 
totics is therefore much slower than in Eq. (f4l|). Since the correction term is negative it is clear 
that the solution to Eq. ( |39|) has a slower growth at low x. That is to say, at low densities solutions 
of both Eq. ( [4l|) and Eq. (|39|) start growing according to the double logarithmic DGLAP equation 

xg oc e ln 1 / x (neglecting the running of a s ). This growth however is tempered faster in the 
solution of Eq. (|3^) than in Eq. ( [i"l|). One could wonder perhaps whether the existence of power 
corrections at large density in Eq. (|39|) is an artifact of the Gaussian truncation of the infinite set 
of equations Eq. ([nj). It is easy to see that this is not the case. The reason is that the mechanism 
of screening in our approach is the appearance of the effective mass proportional to the intensity 
of the chromoelectric field in Eq. (|20|). This type of screening with necessity gives power like 
preasymptotic corrections. It would be very interesting to understand better the relation between 
our present approach and the approach of 

Finally we note that in a recent very interesting paper J23| Makhlin and Surdutovich also 
considered the appearance of effective mass in the context of nuclear collisions. Their approach 
physically is complementary to ours. They consider the emission in the final state in a nuclear 
collision. In this situation the large transverse momentum modes are emitted first whereas the 
low transverse momenta are then emitted and propagate not in the vacuum but in the background 
of these hard modes. The propagation of soft modes in the preexisting hard mode background 
is then characterized by an effective mass. We on the other hand consider fast, large transverse 
momentum fluctuations in the hadron which develop on the background of a slower, smaller 
transverse momentum fields. At large fields the fluctuations are inhibited and this is the origin 
of the effective mass in our approach. It characterizes the emission of hard fluctuations in the 
preexisting soft field background. 
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12 The normalization of the gluon density in [jl^] is slightly different than the one we use throughout this paper. 
This is the origin of the seemingly different factors of 7r and N in Eq. ( fll|) and Eq. (p9|). In fact rescaling g in 
Eq. ( [if] ) by ^ST^ZZu brings the two equations into correspondence in the sense that both, the high density limit 
and the expansion parameter at low density become the same in the two expressions. 
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